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Light-front dynamics can only become a viable alterna- 
tive to the covariant approach if doubts about its covariance 
can be taken away. As a minimal requirement we take that 
the physical quantities calculated with light-front perturba- 
tion theory are the same as those obtained using covariant 
perturbation theory. If this situation occurs, we use the 
word equivalent to characterize it. For quantities that involve 
the calculation of superficially convergent diagrams, proofs of 
equivalence exist. For some types of divergent diagrams the 
proof of equivalence is complicated. Here we deal with dia- 
grams with transverse divergences. Our method is based on 
minus regularization, which is inspired on BPHZ regulariza- 
tion. In a calculation using numerical methods we show how 
to obtain a rotationally invariant amplitude for two triangle 
diagrams contributing to the decay of a scalar boson in the 
Yukawa model. It concludes our proof of equivalence of co- 
variant and light-front perturbation theory. 
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I. INTRODUCTION 

Covariant field theory is the formalism of choice to 
describe situations where creation and annihilation of 
particles are important and where typical velocities are 
comparable to the velocity of light. If the interac- 
tions arc sufficiently weak, perturbation theory is usu- 
ally applied and gives in many cases extremely accurate 
answers. However, in the case of strong interactions, 
or when bound states are considered, nonperturbative 
methods must be developed. Light-front quantization |l]] 
is a Hamiltonian method in which the light-like variable 
x + = (x° + x 3 ) I \/2 plays the role of time, and is therefore 
referred to as light-front time. This method has found 
many applications since it was conceived. Still, some 
problems of a fundamental nature remained. One that we 
are particularly interested in is the question of whether 
full covariance can be maintained in the Hamiltonian for- 
mulation, which is of course not manifestly covariant. A 
partial answer can be obtained in perturbation theory. 
Then the problem can be reformulated as follows: can 
one prove that light-front perturbation theory produces 
the same values of the S-matrix elements as covariant 
perturbation theory? If the answer to this question is 
affirmative, then we use the word equivalent to describe 
the situation. 

The present paper is concerned with one aspect of this 



problem, viz the treatment of transverse divergences in a 
simple model: the Yukawa model with spin-1/2 fermions, 
spin-0 bosons and a scalar coupling. 



A. k -integration and equivalence 

In the work we did before, we used the method of 
Kogut and Soper Q to define light-front perturbation 
theory. This method defines light-front time-ordered 
(x + -ordered) amplitudes by integration of the integrand 
of a covariant diagram, say 



F(q) 



d 4 k I(q;k), 



(1) 



over the light-front energy variable k~ — (k° — k 3 )/^/2. 
In this paper, q always denotes the external momenta 
and k the loop momentum. We can also write Eq. (|l|) 
using light-front coordinates: 

F(q) = J dk+A 2 k^ J Ak-I{q-q\q^-k-k\k^). (2) 

Next, one expresses the integral over k~ , using Cauchy's 
formula, as a sum of residues. One arrives in this way 
at an expression that can be interpreted, possibly after 
recombination of the terms in this sum, as the splitting of 
the covariant amplitude F(q) into a sum of noncovariant 
but light-front time-ordered amplitudes. 

This procedure, sometimes called naive light-cone 
quantization, has been in principle known since the early 
work of Kogut and Soper 01. For convergent diagrams, 
it is nicely pictured in Fig. 111. 
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FIG. 1. The "ideal" case: Outline of our proof of equiv- 
alence of light-front (LF) and covariant perturbation theory 
(PT) for convergent diagrams. The dashed box indicates an 
ill-defined object. 
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The covariant diagram in Fig. [l] is an ill-defined object 
and needs some prescription to give it a definite meaning. 
For example, the measure of the Minkowskian integration 
is not positive definite. The covariant prescription in- 
volves introduction of Feynman parameters to complete 
the squares in the denominator, the removal of terms odd 
in the loop momentum k and the Wick rotation to obtain 
a Euclidian integral. 

It has been the work of Ligterink and Bakker || that 
proves in detail that the rules for constructing light- 
front time-ordered diagrams, explained in many articles 
are correct upon using the /c^-integration prescrip- 
tion. They were the first authors to give a systematic 
derivation of all the different time-ordered diagrams cor- 
responding to a given covariant amplitude, for any num- 
ber of particles involved. If the fc~-integral is convergent 
and the corresponding covariant diagram is also superfi- 
cially convergent, then what remains can be written in 
terms of well-defined, convergent Euclidian integrals. 

When the k~ -integration is divergent, the prescription 
must be altered. Naive light-front quantization fails in 
this case and one must first find a way to regulate the 
/c~-integrals. We proposed in a previous paper || a reg- 
ularization that maintains covariance. There we showed 
that the longitudinal divergences give rise to so-called 
forced instantaneous loops (FILs) and we showed how to 
deal with them such that covariance is maintained. This 
method was also applied to the Yukawa model containing 
spin-1/2 and spin-0 particles. We were able to regularize 
the ^"-integrals for the diagrams with one loop. How- 
ever, in order to show full equivalence to the covariant 
calculation one needs to compute the full integral includ- 
ing the integrations over k + and k . 



B. Ultraviolet and transverse divergences 

Even after the usual procedure has been followed, the 
covariant integral can still be ultraviolet divergent. Lig- 
terink and Bakker did not only discuss diagrams that 
are superficially convergent, but also what to do in cases 
where the covariant diagram is divergent. Their method 
of regularizing divergent diagrams, minus regularization 
is also used in the present paper. A scheme for the 
equivalence of ultraviolet divergent diagrams is given in 

Fig. |. 
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FIG. 2. Outline of our proof of equivalence for diagrams 
with ultraviolet divergences. Dashed boxes indicate ambigu- 
ously defined objects. 

Several techniques are available to remove the ultravi- 
olet divergences, not involving the k~- integration. They 
remain in the light-front time-ordered diagrams as di- 
vergences of the integrals over the transverse momenta. 
Therefore these diagrams are also ill-defined, as indicated 
by the dashed box in Fig. |^. A problem is that many 
of the techniques which are used to regularize covariant 
diagrams have limited use for light-front time-ordered di- 
agrams. For example, one cannot use dimensional reg- 
ularization for the longitudinal divergences. Still, it is 
common to apply it to the transverse divergences. The 
strength of the regularization scheme we use, minus regu- 
larization, is that it does not discriminate between trans- 
verse and longitudinal divergences. Minus regularization 
is based on the Bogoliubov-Parasiuk-Hepp-Zimmerman 
(BPHZ) method of regularization In their pa- 

per, Ligterink and Bakker applied minus regularization 
to three self-energy diagrams. Our contribution is to ex- 
tend their method to more complicated diagrams and 
prove that there is a one-to-one relation between minus 
and BPHZ regularization, such that the physical observ- 
ables found using light-front perturbation theory exactly 
match those found in covariant perturbation theory. 

In the Yukawa model there are five covariant diagrams 
with ultraviolet divergences. The boson and the fcrmion 
self-energy were discussed in our previous article on lon- 
gitudinal divergences. Minus regularization was applied 
and simultaneously removed the longitudinal and the 
transverse divergences. Equivalence was established. 

In two cases we were not able to either find an an- 
swer in the literature or produce ourselves full analytic 
results for the integrals involved; so we had to resort 
to numerical integration. In this paper we discuss these 
two diagrams: the one-boson exchange correction to the 
boson-fermion-fermion vertex and the fermion loop with 
three external boson lines. The first one was considered 
by Burkardt and Langnau E ], who stated that naive 
light-cone quantization leads to a violation of rotational 
invariance of the corresponding S-matrix elements and 
found that invariant results can be obtained using non- 
covariant counterterms. Here we show that no violation 
of rotational invariance occurs if our method of regu- 
larization is applied. Furthermore, our results for the 
light-front time-ordered diagrams sum up to the covari- 
ant amplitude, calculated using conventional methods. 



C. Light- front structure functions 

The two triangle diagrams can be written in the form 
of a sum of tensors in the external momenta, multiplied 
by scalar functions, which we call (covariant) structure 
functions. After splitting a covariant diagram in light- 
front time-ordered ones, these can be written again in 
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terms of tensors multiplied by functions of the external 
momenta. The latter are called light-front structure func- 
tions. They are not invariant as they are not defined by 
four-dimensional invariant integrals, but rather by three- 
dimensional integrals. The different structure functions 
have different divergences and they must be treated ac- 
cording to their types of divergence, which we enumerate. 

1. Light-front structure functions without transverse 
divergences. Neither the covariant nor the light- 
front formulation contains any divergences. Inte- 
gration over k" suffices to prove equivalence. Mi- 
nus regularization is not allowed. 

2. Light-front structure functions with cancelling 
transverse divergences. The individual light- 
front time-ordered diagrams contain divergences 
not present in the covariant amplitude. Applica- 
tion of minus regularization to the time-ordered di- 
agrams is not allowed. We show that the diver- 
gences cancel if all the time-ordered diagrams are 
added, and that their sum equals the corresponding 
covariant amplitude. 

3. Light-front structure functions with overall trans- 
verse divergences. Divergences appear in the co- 
variant amplitude as well as in the light-front time- 
ordered diagrams. We apply BPHZ regularization 
to the covariant amplitude and minus regulariza- 
tion to the time-ordered diagrams. 

For the first two cases one can prove equivalence using 
analytic methods alone. The proof of equivalence can 
be found in Refs. H||. For the structure functions with 
overall transverse divergences we have to use numerical 
techniques. We show that for the decay of a boson at 
rest, for both triangle diagrams, one obtains a rotational 
invariant amplitude, identical to the covariant calcula- 
tion using BPHZ regulariztion. The fifth diagram with 
transverse divergences, the fermion box, will not be dis- 
cussed. 



D. Outline 

The setup for this article is as follows. In Sec. [il] we 
introduce minus regularization. In Sec. [II and Sec. |y| 
we discuss the equivalence of covariant and light-front 
perturbation theory for the fermion triangle and the 
one-boson exchange correction. In both cases we start 
with the covariant calculation and do the BPHZ regu- 
larization if necessary. Then we calculate the light-front 
time-ordered diagrams and apply the method mentioned 
above. In both cases, we conclude by giving a numerical 
example of rotational invariance. 



II. MINUS REGULARIZATION 

Minus regularization is inspired by the BPHZ method 
of regularization, which gives finite and covariant results. 
By construction, we ensure that minus regularization 
does the same. First we sketch the method in the case 
of one-loop diagrams with one independent external mo- 
mentum (self-energies), and next when two independent 
external momenta (triangle diagrams) are present. We 
conclude by generalizing this to a one-loop diagram with 
n external momenta. For convenience, we shall assume 
in the latter case that only logarithmic and linear diver- 
gences are present, such that only the first term of the 
Taylor expansion around the renormalization point needs 
to be subtracted. 

Wherever we use the word "amplitude" in this section, 
we refer to an invariant function of the external momenta. 
It is understood that the integrals defining the invari- 
ant functions are formally written down in terms of four- 
dimensional integrals, which are split into time-ordered 
pieces by integration over k~ . 



A. One external momentum 

First we discuss the simple case of one external mo- 
mentum, which can be applied for self-energy diagrams. 



1. BPHZ regularization 

We start with the BPHZ regularization method, which 
can be applied to covariant diagrams. The amplitude has 
the following form 



d 4 fc J cov (^;A;) 



F(0) + q 2 F'(0) 



(3) 



where I cov (q 2 , k) is the covariant integrand generated by 
applying standard Feynman rules. BPHZ regularization 
renders the amplitude finite by subtracting the infinite 
parts. We choose the point q 2 = as the renormalization 
point, around which we expand the amplitude in a Taylor 
series. The higher orders in the expansion (^) are denoted 
by the elipsis. The regularized amplitude is then 



F R (q 2 ) = F(q 2 ) - F(0). 



(4) 



However, this is a purely formal operation, since we are 
subtracting two infinite quantities. It is better to write 



(5) 



F R (q 2 ) = J d 4 k (L cov (q 2 ; k) - I cov (0; k)) 

= f dq> 2 Jd*k^i cov ( q ' 2 ;k). 



This guarantees that the amplitude becomes finite. 
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2. Minus regularization 

Typical for minus regularization is that one writes the 
amplitude, as well as the renormalization point, in light- 
front coordinates. The covariant choice q 2 — corre- 
sponds to q~ = q /{2q + ). A time-ordered amplitude 
corresponding to the covariant form (|^) can be written 
in light-front coordinates as follows 



F(q- q+ q X ) = J d 3 k I lho (q- q+ q X ; k) 
F(^, q + , q X ) + 2q+(q-- ^_)F'(^, q+ q X ) + 



(G) 



2q+' K 2q^ 



where /ift is the integrand of the light-front time-ordered 
diagram, which was generated by integrating the covari- 
ant integrand J cov over k~ as is explained in Ref. ^). The 
prime denotes differentiation with respect to q~ . Similar 
to Eq. (|^) we can write the regularized amplitude as 



F R (qlql qi -q 2 ) = J d 4 fc 

x (lc OV (qj,ql,qi-q2;k) - l cov (0;k)) 



(10) 



We cannot, as in the previous section, differentiate with 
respect to all external momenta. We would then subtract 
finite parts from the Taylor series, containing physical 
information. This can be circumvented by introducing a 
dummy variable A, which parametrizes a straight line in 
the space of the invariants between the actual external 
momenta q\,q\,q\ - q2 and the renormalization point: 



F R (qlqlqi-q 2 )= J dA J d 4 fc 
d 

X ^/ C ov(Agi, Xq 2 , Xqi-q 2 ; k) 



(11) 



We have verified that the A-method gives the correct re- 
sult for the case where one independent external momen- 
tum occurs. 



F MR (q-q+q^ = J^dq'-J d 3 fc -^-I^q'-q+q^k). 

2q+ 

(7) 

So far we have described the minus regularization method 
introduced by Ligterink and Bakker B . 



B. Two external momenta 

In Ref. three self-energy diagrams were discussed. 
For the triangle diagram the minus regularization tech- 
nique needs to be extendec^]. We will tune the technique 
by comparing it to BHPZ regularization. 



1. BPHZ regularization 
The amplitude has the following covariant form 

F{ql,ql,qx-qi) = J d 4 fc I cav (qf, ql, qi-qr, k) 
= F(6) + qfF{(0) + q 2 2 F^6) + qi -q 2 F^0) + ... (8) 

where is the renormalization point q\ = q\ = q±-q 2 = 
and F[ is the derivative of F with respect to the ith 
argument. 

F R (qlql qi -q 2 ) = F(<£, ql q v q 2 ) - F(0). (9) 

Again, this is a purely formal operation, since we are 
subtracting two infinite quantities. We write 



2. Minus regularization 

Again, we write the amplitude in the light-front time- 
ordered case as a three-dimensional integral: 

F(q~, qf, qt) = J ^ k ^(q^ qf, ql\ k )- (12) 
The regularized amplitude is 

F R (q-, q+ qi) = F(q7, q+ - F{r-, r+ r±), (13) 

where r defines the renormalization surface. It is a hy- 
persurface determined by the following conditions: 



r? 


= 2r^r\ 


-^ 2 = o, 


rl 




- 4 2 = o, 


T\-r 2 


= r i r i h 





This set of equations is equivalent to 
r\ = 0, r 2 = xn. 



(14) 



(15) 



The rf enter in the integration boundaries; therefore we 
would like them to remain unaffected by regularization 
(jl = ql). This implies that x can be found from 



4 

x = — ■ 



(16) 



The only freedom that remains is the choice for r ^ . Two 
choices come easily to mind: rf- = (method MRO) and 
= q^ (method MR1). 



1 We suggest the name MR + . 



(MRO) 
(MR1) 



ri = 0- 1 => ri = X , (17) 
rf = qi => r x = X qi- (18) 
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TABLE I. The light-front parametrization of the renor- 
malization point r M for two equivalent choices of minus 
regularization, MRO and MR1. 





MRO 


MR1 




(0,9+0^) 






X(0,g+ 0^) 





The details are worked out in Table I. The light-front 
coordinates of the renormalization point are used in the 
following way to find the regularized light-front ampli- 
tude: 

F^(qr, q +, q t)=j\xjd 3 k^ 

x /ifto(A(«T- U)+rr,q+ \(qj- rf)+rf ; fc). (19) 
In this formula we recognize our choice rf — qf . 



regularization, the renormalization point corresponds to 
all invariants connected to the external momenta being 
equal to zero. These conditions allow minus regulariza- 
tion to take on a number of forms. Of these, we shall ap- 
ply MRO and MR1. The main difference between them 
is that MRO does not choose one of the momenta as a 
preferred direction, and therefore it explicitly maintains 
all symmetries of the external momenta. Furthermore, 
MRO gives rise to shorter formulas for the regularized 
integrands. 

In the next two sections both methods are being ap- 
plied to the parts of two light-front time-ordered trian- 
gle diagrams in the Yukawa model containing transverse 
divergences, viz the fermion triangle and the one-boson 
exchange correction. 

III. EQUIVALENCE FOR THE FERMION 
TRIANGLE 



C. Several external momenta 

The method just described can be generalized to the 
case of a loop with an arbitrary number of external lines. 
The procedure is almost the same as for two external 
momenta. The renormalization surface is given by 



rf = 2r7r+ - 2 = 0, 
■rj =rr r + +r+rj - r; t 



These equations are equivalent to 



(i^j). 



rf=0, 



X^i- 



(20) 
(21) 



(22) 



Again, we make the choice to leave the plus components 
of the momenta unaffected by regularization: rf = qf. 
This implies that the x% are fractional longitudinal light- 
front momenta. 



2L 
it 



(23) 



Two choices for are listed below. This then determines 
all other rr- 



(MRO) 
(MR1) 



= J 



= J 



1i Ti = XiQi 



(24) 
(25) 



D. Summary 

The way we setup minus regularization does not rely 
on the structure of the covariant or the time-ordered 
diagrams, but works on the level of the external mo- 
menta only. If an amplitude has a covariant structure be- 
fore regularization, minus regularization guarantees that 
it remains covariant. In our implementation of BPHZ 



In the Yukawa model there is an effective three boson 
interaction, because for a fermion loop with a scalar cou- 
pling Furry's theorem does not apply. The leading order 
contribution to this process is the fermion triangle. A 
scalar boson of mass fi and momentum p comes in and 
decays into two bosons of momentum q\ and q2 respec- 
tively. The fermions in the triangle have mass to. The 
covariant expression for the amplitude is 



d 4 fc Tr Pi + m)(#2 + m)( j£ + m)} 



Mi: 



(k\ — to 2 )(/c| — m 2 )(k 2 — to 2 ) 



(26) 



The subscript "Min" denotes that the integration is over 
Minkowski space. The usual imaginary parts of the Feyn- 
man propagators have been dropped. We have omitted 
numerical factors and have set the coupling constant to 
unity. The momenta k\ and k 2 indicated in the diagram 
are given by 

k\ = k-q\, k 2 = k + q 2 . (27) 
Of course, by momentum conservation we have 



p = qi + <?2- 



(28) 



We evaluate the integral (|26|) first in the usual covariant 
way, and subsequently carry out ^"-integration to pro- 
duce the light-front time-ordered diagrams. Note that 
integral (|26|) is an ill-defined formula. In both methods 
mentioned we have to define what we mean by this inte- 
gral. 
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A. Covariant calculation 

The following method is usually applied to calculate 
the fermion triangle in a covariant way. First, one in- 
troduces Feynman parameters x\ and x%, and then one 
shifts the loop variable k to complete the squares in the 
denominator. The result is 



1 r~L—xx 
dx\ I dx 2 
(i Jo 



d 4 k 



(29) 



3fc 2 



Min 

V 2 ) + terms odd in k 



(k 2 



Q 2 



with 



Q 2 = Xi(l-xi) q\ + x 2 (l-x 2 ) q\ + 2xix 2 qrQi, (30) 
V 2 = xi (3xi- 2)q 2 + x 2 (3x 2 -2)q 2 

+ (2(xi+x 2 ) - 6xix 2 - 1) q%-q 2 . (31) 

As a last step, we remove the terms odd in k. 

B. BPHZ regularization 

The regularized fermion triangle can be found by ap- 
plying the BPHZ regularization scheme ([ll]) to the co- 
variant formula (^). The integral is now finite; so we 
can do the Wick rotation and perform the k integrations. 
The result is 











i = — 4ir 2 i / d:ci 


f dx 2 [ 


JO o 


'o Jo 



m(m 2 (5Q 2 -V 2 



dX 
6AQ 4 ) 



(to 2 - \Q 2 ) 



(32) 



The superscript R indicates an integral regularized ac- 
cording to the BPHZ method. 



C. Light- front calculation 

Using the method given in Ref. Q we proceed as fol- 
lows. The k~ dependence of a spin projection in the 
numerator is removed by separating it into an on-shell 
spin projection and an instantaneous part: 



where the vector /cf on is given by 

.+ u±\ - ( k ^ + m2 u+ i,A 



[ k i j ov: 



2W 



(33) 



(34) 



Factors like (fc — fc i on ) can be divided out against prop- 
agators and this cancellation gives rise to instantaneous 



fermions. The integration over k is performed by con- 
tour integration. The poles of the propagators are given 
by 



H~ = 



H 



2k+ 



H 2 =-<h 



ki + rn 2 
2k+ ' 

k^ 2 + m 2 

-| £ . 



(35) 
(36) 

(37) 



This integration gives rise to the different time-ordered 
diagrams, as explained in more detail in Refs. |,|. The 
result is 





(38) 



The diagrams on the right-hand side are light-front time- 
ordered diagrams. Time goes from left to right. The pic- 
tures can be recognized as time-ordered diagrams because 
of the time-ordering of the vertices and the occurrence 
of instantaneous fermions, indicated by a horizontal tag. 
Explicitly; 



2m J d 2 k^ 
Trp lon + 



it 



dfc H 



o 8kfk^k+ 

TO) (^ 2 Oil O 



(H^-H 2 )(H^-H-) 



(39) 




= 2m / d 2 k 



dk+ 



, a 8k+k+k+ 

Trp 

Ion on 

= 2m / d 2 k^ 



(40) 



dfc H 



Trp 

Ion + rn) ( # 2 on 

+ m)7 + ] 

H 1 -H 2 



(41) 



= -2m / d 2 k 



dk+ 



. q + 8k+k+k+ 
Tr[(fr 

Oil 

+ TO) (^2 on on + to)] 



(42) 



G 




2ni / cYk 



2 I.JL 



Tr[ 7 +(# 

2 on + m)()t o 1.1 + m)] 
H~-H~ 



(43) 



2ni / d 2 fc J 



_ g + 8k+k+k+ 



Tr[(& 

on + m) ( ^2 on 

+ to) 7 j 

■Hi — -ffo 



(44) 



Note that the diagrams ( |4l| ) and (EJ) with the instan- 
taneous exchanged fermions have the same integrand. 
However, the longitudinal momentum k + has a different 
sign. 

Although we could have expected diagrams with two 
instantaneous fermions, we see that they are not present. 
This is so because we use a scalar coupling and there- 
fore two 7 + matrices becoming neighbors give 0. No 
so-called forced instantaneous loops are present. These 
FILs obscure the equivalence of light-front and covariant 
perturbation theory and have been analyzed in Ref. 
They will not be discussed in this paper, since they are 
related to longitudinal divergences. 

The traces can be calculated. We obtain 



Ti-p lon +m)(# + m)(ji 

on + to)] 

= 4to(to 2 + k lon -k on + k 2on -k on + k lon -k 2on ), (45) 
Tr[(# lon + m)(# 2on + to) 7 +] = 4to (2fc+- q++ q+), (46) 
Tr[(# lon +TO) 7 +(# on +TO)] =4m(2k+-q+), (47) 
Tr[ 7 +(#2on +m)(# on +to)] =4TO(2fc++ g +). (48) 

We see that the high orders in k 1 - have disappeared in the 
traces. However, logarithmic divergences remain in all 
light-front time-ordered diagrams (p9|)-(^4]). We tackle 
them with minus regularization, as introduced in the pre- 
vious subsection. 



D. Equivalence 

As the fermion triangle is a scalar amplitude, there is 
only one structure function present. It belongs to the first 
category we mentioned in the Introduction: it is logarith- 
mically divergent, but has no longitudinal divergences. 



1. Light- front structure functions with transverse divergences 

We applied minus regularization to the integrands of 
the six light-front time-ordered diagrams, using both the 
MRO and MR1 methods. We used mathematica to 
do the substitution and the differentiation with respect 



to A, given by Eq. (|l9|). However, mathematica was 
not able to do the integration, neither analytically nor 
numerically. Therefore the integrand was implemented 
in FORTRAN which was well capable of doing the four- 
dimensional integration using imsl routines based on 
Gaussian integration. 

Because the integrations cannot be done exactly, we 
saw no possibility of giving a rigorous proof of the equiv- 
alence of light-front and covariant perturbation theory. 
Instead we make a choice for the parameters, such as 
the masses and the external momenta, and show that 
our method gives the same result as the covariant cal- 
culation with BPHZ regularization. We calculated the 
decay amplitude of a scalar boson at rest, as is pictured 
in Fig. | 




FIG. 3. A boson is at rest and decays into two particles 
flying off in opposite directions. The angle 8 is the angle 
between the momentum of one of the fermions and the z-axis. 

From a physical point of view, there is no preferred 
direction, and therefore we demand that our choice of 
the coordinates of the light-front have no influence on 
the outcome of the calculation. The decay amplitude, 
which is a scalar quantity, should give the same result 
for each possible direction in which the bosons can fly 
off. 

There are six minus-regularized light-front time- 
ordered fermion triangle diagrams contributing to the 
boson decay. Each individual light-front time-ordered 
diagram has a manifest rotational invariance in the x- 
y-plane, and therefore we expect the same for the sum. 
However, since light-front perturbation theory discrimi- 
nates between the z-direction and the other space-like di- 
rections, the light-front time-ordered diagrams can (and 
should) differ as a function of the angle, 9, between the 
momentum of one of the particles flying off and the z- 
axis. The absolute value of the momentum was fixed. It 
is not immediately clear that the sum should be invari- 
ant. This investigation becomes more interesting since 
it is believed |lj] that rotational invariance is broken 
in naive light-cone quantization of the Yukawa model. 
However, the results shown in Figs. demonstrate that 
rotational invariance is not broken. Note that we have 
dropped the factor — i common to all diagrams. 



7 




FIG. 4. The thick line at a value of 125 represents the sum 
of the six light-front time-ordered amplitudes. It is indepen- 
dent of the angle 9, defined in the previous figure. The four 
largest contributions come from the diagrams without instan- 
taneous parts (solid lines) and the diagrams with an instanta- 
neous exchanged fermion (dashed lines), as indicated by the 
diagrams. 

Two light-front time-ordered diagrams (|4C|) , (Jl3|) con- 
tributing to the boson decay and indicated by double- 
dashed lines are so small they can hardly be identified in 
Fig. ||. In Fig. || we depict these two on a scale that is 
a factor of 100 larger. In the same figure we show the 
difference of the sum of the six light-front time-ordered 
diagrams (using MR1 and 128 points in every integration 
variable) and the covariant result. It has a maximum of 
0.03%. 



-0.5 



/ 



\ 



tt/4 it/2 3tt/4 6 -> it 

FIG. 5. The amplitudes of the two small contributions 
(double-dashed lines) and the difference between the sum of 
the six light-front time-ordered diagrams and the covariant 
amplitude (thick solid line). 

In Figs. ^ and ^| we see that interchanging the outgoing 
bosons is the same as replacing 8 by tt — 8. 

We verified that the individual diagrams are rotational 
invariant around the z-axis. We illustrate this in Fig o. 



MR 





MR 




(xlOO) 



+ 



MR 



+ 



BPHZ 




FIG. 6. Commutative diagram of the boson decay ampli- 
tude. The boson is at rest in the origin and decays. The 
outgoing bosons fly off in opposite directions. Points on the 
surfaces have polar coordinates (A, 9, <f>), where A the is mag- 
nitude of the amplitude and 9 and <f> are the polar angles of 
the momentum of one of the outgoing particles, as defined 
in Fig. |^. Because the diagrams on the second line are very 
small, the scale has been enlarged by a factor of 100. For 
the light-front time-ordered diagrams on the first three lines 
minus regularization (both MR0 and MR1) is used, for the 
covariant diagram on the last line we used BPHZ regulariza- 
tion. 

Summing up, we find that the sum of the minus reg- 
ularized light-front time-ordered diagrams is rotational 
invariant. The deviation from the covariant result is 
smaller than 0.03%. It is illustrated in Fig. [|. We 
checked, by varying the number of integration points, 
that the deviations are due to numerical inaccuracies 
only. We conclude that, for the fermion triangle, the 
covariant calculation in combination with the BPHZ reg- 
ularization scheme gives the same result as the light-front 
calculation in combination with minus regularization. 
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IV. EQUIVALENCE FOR THE ONE-BOSON 
EXCHANGE CORRECTION 

The second process under investigation was studied 
before by Burkardt and Langnau p2 [. A scalar boson 
of mass yU and momentum p decays into two fermions 
of mass m and momentum q\ and q 2 respectively. The 
lowest order correction to this process is the one-boson 
exchange correction. The amplitude is given by the inte- 
gral 




d 4 fc (#1 + m)(#2 + m) 
(kf — m 2 )(k 2 — m 2 )(fc 2 — /i 2 ) 



(49) 



Again, this equation is undefined as it stands. First we 
have to make it a well-defined object. In Sec. IV A we 



apply the covariant method and in Sec. IV C we use light- 
front coordinates. 



(56) 



The functions F l depend on the masses and the exter- 
nal momenta q 2 , q\ and q\-q 2 . If we define the integral 
operator 

/[/] = 2 / da* / ^ dx 2 f d 4 fc (fc 2 - M 2 + Q 2 )" 3 /, 

JO JO ./Min 

(57) 

then we have, using 4i = q 2 , etc. 



F 1 = I [fc 2 + m 2 - V 2 ] , 
F 2a = 2ml [1-xi- x 2 ] , 
F 2s = 2ml [-xi + x 2 ] , 

F 3 = I[l-x 1 -x 2 ]. 



(58) 
(59) 
(60) 
(61) 



We see that the only function which needs to be regular- 
ized is F 1 . The functions F 2 and F 3 are convergent and 
do not require regularization in a covariant calculation. 



A. Covariant calculation 



B. BPHZ regularization 



Using Feynman parametrization the one-boson ex- 
change correction can be rewritten as 



pi fl— Xl p 

= 2 1 dxi / dx 2 / d 4 fc 

Jo Jo JMin 



(50) 



ft 2 + ((1-aci) di+x 2 ) 42 + m) (-xi 4i-{l-x 2 ) 4 2 + m) 



(fc 2 - M 2 + Q 2 



terms odd in fc 



with 



M 2 = (xi + x 2 )m 2 + (1 - xi - x 2 )^ 2 , (51) 
Q 2 =xi(l- xi)qf + x 2 (l - x 2 )q 2 + 2x\x 2 q v q 2 , (52) 

and where terms odd in fc in the numerator are not speci- 
fied, since they will be removed according to the covariant 
prescription. We also define 



V 2 = Qf + {l-xi-x 2 )q v q 2 . 



(53) 



^From Eq. (Q) we can infer that the Dirac structure 
of the diagram is 



= F 1 +F 2 ^+F^\[ 1 ^ lv \. (54) 



where the vector part contains a symmetric and an anti- 
symmetric part, 

F 2 » = F 2s (tf + <g) + F 2 *(q» - q»l (55) 

and the tensor part has the form 



The regularized structure function F 1R can be found 
by applying the BPHZ regularization scheme (|ll]) to the 
structure function (|5|). The integral is now finite; so we 
can do the Wick rotation and perform the fc integrations. 



F 1R (q 2 ,q 2 , qi -q 2 ) = -2ir 2 i Ux t f dx 2 f dX 

Jo Jo Jo 

(Q 2 (XP 2 -m 2 ) 



Q 2 + \V 2 



\2{M 2 -\Q 2 ) 2 M 2 -XQ 2 



(62) 



We have not been able to do all three integrations exactly. 
The A integration and one of the x integrations can be 
done analytically, and the remaining integration numer- 
ically. As F 2 ^ and F 3 do not need to be regularized, 
this concludes the covariant calculation of the one-boson 
exchange correction. 



C. Light- front calculation 

In our previous paper [|| it was shown how to derive 
the light-front time-ordered diagrams corresponding to 
the covariant diagram (|4^) using fc _ -integration. One can 
write the time-ordered diagrams individually, or one can 
combine propagating and instantaneous parts into so- 
called blinks. Blinks, introduced by Ligterink and Bakker 
||, have the advantage that the l/fc + -singularities cancel 
and the number of diagrams is reduced. 

In the two triangle diagrams studied here it makes no 
difference whether blinks are used or not. In the case 
of the fermion triangle we calculated light-front time- 
ordered diagrams. Here we use blinks, to demonstrate 
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that our technique also works in this case. The one-boson 
exchange correction has two blinks: 



(63) 



The poles of the two fermion propagators in the tri- 
angle are given by Eqs. ([36]) and (|37]). The pole of the 
boson propagator is given by 




H' 



k 



2k+ 



(64) 



The amplitudes including blinks are 



"\ = -2m / d 2 fc 



dfc H 



on 

+ m) 

(h--h;)(h--h;) ' 



(65) 



2ni I d 2 ^ 



io Skjkjk+ 
(faon + m)(}i lon + tf + m) 

We will now focus on the blink in Eq. (|66|). 

because we can use 

Therefore we obtain 



(66) 

It simplifies 
(67) 



2m I d 2 ^ 



«i dk+ 



o Sk{kjk+ 



2m 2 + )i lon (0 + 2m) 
(H--H-)(H--H-) 



(68) 



In the same way as we did for the covariant amplitude 
we can identify the different Dirac structures 



Fl+F^+F^l^, 7 „}. (69) 



Although at first sight it looks as if the diagram in 
Eq. (|6^) has a covariant structure, covariance is spoiled 
by the integration boundaries for k + . Therefore these 
functions are not covariant objects. We have to investi- 
gate equivalence for the structure functions separately. 

The light-front structure function F± can be found by 
taking the trace of Eq. (|6^) , since all the other structures 
are traceless. Carrying out the traces one finds 



F{ = 2ni / d 2 k 



^ dk+ 



2m 2 + k lon -p 

The other structures of the blink diagram ( p8| ) are 

<-*t dk+ 



(70) 



F 2m - 2ni ( d 2 fc^ ( 1 — 
1 - 2m J dk J 8 ki 



8fc 1 + fc+/^^ 



2m (k loa Y 



(H^-Hz)(H^-H- 



(71) 



1 J Jo 8k+k+k+ 



(fclon) M V V 



(H£-H^)(Hr-H-) 



(72) 



In a similar way we can derive the structure functions 
corresponding to the other blink diagram. 



D. Equivalence 

We can identify the different types of divergences, as 
explained in the Introduction. 



1. Light-front structure functions without transverse 
divergences 

The parts of the blinks without any ultraviolet diver- 
gences are F 2f * and F^" , except for \i being — . No 
cancellations need to be found and no regularization is 
necessary. 



2. Light-front structure functions with cancelling transverse 
divergences 

In the last two structure functions we see something 
odd happening. Both F 2 ^ and F i >J ' v are divergent for /i 
being — . However, these divergences are not present in 
the covariant structure functions F 2tL and F 3mi/ . It would 
we illegal to apply minus regularization, since the covari- 
ant amplitude does not need to be regularized. We found 
that the divergences corresponding to the first blink can- 
cel exactly against those of the second blink. To simplify 
the calculation we use internal variables x' and k and 
external variables x, 0.7 an d it ■ These are introduced in 
Appendix |A|. 

We have to verify the following relation of equivalence 



F 2 - = F 2 - 



Fo 



(73) 



According to the reasons mentioned above we have to 
demand that the divergent parts in the right-hand side 
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cancel. We find that only the highest order contribution 
in k contributes to a divergent integral, because we can 
write 

= j&^f^p+ft^))* (74) 

where gf~(k ± ) is the part of the integrand without ul- 
traviolet divergences, and the term with ft gives rise to 
a logarithmically divergent integral. We have to check if 
we have 

ft + ft = 0- (75) 

In Appendix ^ the full formulas for the functions ft are 
given, from which it follows that condition ( j75| ) holds. 
For /i being — in the structure function Ff^ one can 
apply the same method. 

3. Light-front structure functions with transverse divergences 

The structure function F 1 in the covariant calcula- 
tion contains an ultraviolet divergence. In the light-front 
structure functions these appear as divergences in the 
transverse direction. The equation under investigation is 
the following: 

F 1MR + F 1MR = F 1R (?6) 

For the same reason as for the fermion triangle, an an- 
alytic proof of this equation is not possible. We investi- 
gated rotational invariance of the left-hand side of this 
equation, and furthermore we checked if it gives the same 
result as the covariant calculation on the right-hand side. 
A boson is at rest and decays into two fermions as indi- 
cated in Fig. ^. The fermion mass is taken to be the same 
as the boson mass. Therefore there can be no on-shell sin- 
gularities of intermediate states. Also, we dropped the 
common factor —i. The contributions of the two blink 
diagrams are given in the commutative diagram of Fig. |t]. 
We made the arbitrary choice of applying minus regular- 
ization MR1, and used 128 points in integration variable. 

The error, i.e., the difference between the covariant cal- 
culation with BPHZ regularization and the sum of minus 
regularized blinks, has a maximum of 0.02%. This devia- 
tion results from numerical inaccuracies, as was checked 
by varying the number of integration points. 

We conclude that no significant deviation from a ro- 
tational invariant amplitude is found. Moreover, we 
found that the sum of the light-front time-ordered di- 
agrams is the same as the covariant amplitude for the 
one-boson exchange correction. Again, the procedure of 
^"-integration and minus regularization proved to be a 
valid method. 



z 




FIG. 7. Commutative diagram of the one-boson exchange 
correction. A boson decays at rest. The outgoing fermions fly 
off in opposite directions. The distance from the origin gives 
the amplitude of the regularized diagram for the fermion fly- 
ing off in this direction. For the light-front structure functions 
on the first two lines, minus regularization (MR1) is used; 
for the covariant structure function on the last line, we used 
BPHZ regularization. 

V. CONCLUSIONS 

In the Yukawa model with a scalar coupling there are 
five single-loop diagrams with transverse divergences, of 
which two also contain longitudinal divergences. For 
all other one-loop diagrams and all multiple-loop dia- 
grams that do not contain subdivergences, the proof of 
the equivalence of covariant and light-front perturbation 
theory was given by Ligterink and Bakker || upon using 
the ^"-integration prescription. For the two single-loop 
diagrams with longitudinal divergences this integration 
is ill-defined. This problem was dealt with in a previous 
paper g. 

Of the three remaining diagrams two are thoroughly 
analyzed in this paper. For the parts of these diagrams 
without transverse divergences the ^"-integration recipe 
of Ligterink and Bakker applies. For the parts with trans- 
verse divergences a proof of equivalence is complicated by 
the fact that the amplitudes depend on three independent 
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scalar products of the external momenta. We applied an 
extended version of the method of minus regularization 
invented by Ligterink and Bakker. It is on a friendly 
footing with the light-front, because it can be applied to 
both longitudinal and transverse divergences. Moreover, 
it has strong similarities to BPHZ regularization, which 
is suitable for covariant perturbation theory. We were 
able to tune the regularization in such a way that mi- 
nus regularization is analogous to BPHZ regularization. 
Therefore, we expect an exact equality between the co- 
variant and the light-front amplitudes. We showed that 
rotational invariance is maintained and we expect that 
other nonmanifest symmetries on the light-front, such as 
boosts in the x-y-plane, are also conserved. 

The final formulas obtained did not yield to analytic 
integration. Therefore we had to resort to multidimen- 
sional numerical integration. As rotational invariance 
was shown previously to be violated in naive light-cone 
quantization fl^| , we investigated rotational invariance, 
which is one of the nonmanifest symmetries on the light- 
front. Our results demonstrate, within the errors due to 
the numerical methods used, that covariant and light- 
front time-ordered perturbation theory give the same 
physical matrix elements. 

One diagram with transverse divergences has not been 
discussed in our two papers on equivalence, namely the 
fcrmion box with four external boson lines. It is a scalar 
object, similar to the fermion triangle. The results ob- 
tained for the latter convinced us that upon minus reg- 
ularization we shall find a covariant result. As there are 
more time-orderings, and because one cannot test for ro- 
tational invariance as easily as for the triangle diagrams, 
we did not investigate this much more complicated situ- 
ation. 

We trust that with our elaborate discussion of diver- 
gent diagrams in the Yukawa model we have illustrated 
the power of minus regularization and taken away doubts 
about the covariance of light-front perturbation theory. 
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k+ 

x = — = (x- l)x, 



k + + q 2 _ x + x 

qt x 



(Al) 
(A2) 



Or, equivalcntly, 



k + = x'ql = {x- l)q£, 
kf = (x' - l)q+, 
k 2 = xq 2 . 

In the numerator of the integrals defining light-front 
structure functions we encounter on-shell spin projec- 
tions. They can be rewritten in terms of internal vari- 
ables using 



^lon ^ 



^2on — 



ki + m 2 
2(x' - l)q+ 

2 

2xqt 



(A3) 
(A4) 



The energy denominators can also be written in terms of 
internal and external variables. The poles are given by 
Eqs. @, © and d§): 

2q+(H 1 -H 2 ) = 2q+ [p +^ - > 



2k1 



2k 2 



= (P 2 +P ±2 ) 



2q+(H--H-) = 2q+ q- 



-i i I 9 i I 9 

j_2l+X »i +m k 2 +m 



1-x' 



(A5) 



k +fx ^i +n / 



2k+ 

„-L z i ,,2 



2kJ 



. I 2 o , I 2 n 

ii+ii zr 1 — — ( A6 ) 



x 



1-x' 



2q+(H--H-) = 2q+ q~ 



k 1 - +i? k^ +m 2 



2k+ 



i j 2 k^+fi 2 k 2 2 +m 2 . 

= + 4 1 • (A7) 

1 — x x 

The integration measures can be rewritten as follows: 

rl Ax' 



2m I ; — r — 7-, = —71"* 



o 8k+k+k+ ' Jo (I- x')xx' ' 

1 da; 



. f° dfc+4q+g+ 



(l-x')x(l-x) 



(A8) 



(A9) 



APPENDIX A: INTERNAL AND EXTERNAL 
VARIABLES 

We get more insight into the properties of the struc- 
ture functions if we rewrite them in terms of internal and 
external variables. This can be done by defining 



We conclude that it is possible to write the struc- 
ture functions in terms of the external variables 
qi , q 2 , q± , <?2~ an d X an d integrals over the internal 
variables x or x' and k 1 - . The divergent part of the struc- 
ture functions F 2 can now be written as 
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ft 



1 



dx' to q[ 

(l-x')xx' (x>-l)q+ ^ 
i/i i \ -1 



1 - X? 

•2- 



1 

xx 



1 



-7TZ 



1 

x 1 — X 
dx 



o (1 — x')x(l — x) xq\ 



m 



1 



1-x' 



1 



1 -x 



(A10) 



(All) 



Upon cancelling common factors, and using Eq. (A2), we 
can evaluate the integrals and obtain 



,2- ,2- ■ X rn .to 
Ji =-/2 =7r *-;— + =7r *-r- 

Therefore condition is verified. 



(A12) 
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